The metallicity of the old open cluster NGC 6791 
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ABSTRACT 

We have observed four red clump stars in the very old and metal-rich open cluster NGC 6791 
to derive its metallicity, using the high resolution spectrograph SARG mounted on the TNG. 
Using a spectrum synthesis technique we obtain an average value of [Fe/H] = +0.47 (±0.04, 
r.m.s. = 0.08) dex. Our method was tested on /i Leo, a well studied metal-rich field giant. 
We also derive average oxygen and carbon abundances for NGC 6791 from synthesis of [O i] at 
6300 A and C 2 at 5086 A, finding [O/Fe] ~ -0.3 and [C/Fe] ~ -0.2. 



Subject headings: Stars: abundances 
individual: NGC 6791 

Introduction 



Galaxy: disk - Galaxy: open clusters - Galaxy: open clusters: 



The determination of the abundances in stars 
of different age and Galactic location is one of the 
basic tools to interpret the chemical evolution of 
the Milky Way disk. Galactic open clusters are 
particularly well suited to this purpose (e.g., Friel 
1995) since they reach ages as old as the disk itself, 
cover a large range in metallicities and ages, are 
observed in regions of the Galactic disk likely char- 
acterized by different star formation histories, and 
their distances and ages can be determined with 
a precision not reachable for field stars, except for 
the nearest ones. 

We are presently studying in an accurate and 
homogeneous way a significant sample of open 
clusters (Bragaglia & Tosi 2006 and references 
therein): reliable distances, reddenings and ages 
are derived from photometry with the synthetic 



1 Based on observations made with the Italian Telesco- 
pio Nazionale Galileo (TNG) operated on the island of La 
Palma by the Fundacion Galileo Galilei of the INAF (Isti- 
tuto Nazionale di Astrofisica) at the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astrofisica 
de Canarias 

2 INAF-Osservatorio Astronomico di Padova, vicolo Os- 
scrvatorio 5, 1-35122 Padova, Italy 

3 INAF-Osservatorio Astronomico di Bologna, via Ran- 
zani 1, 1-40127 Bologna, Italy 



color-magnitude diagrams technique (Tosi et al. 
1991), and will be combined with metal abun- 
dances from high resolution spectroscopy to pro- 
vide robust constraints on the current and past 
disk properties. 

We have already presented the detailed chem- 
ical abundances of a few old open clusters (Bra- 
gaglia et al. 2001; Carretta et al. 2004, 2005). To 
our sample we now add NGC 6791 which, with an 
age of about 9-10 Gyr, is one of the oldest open 
clusters of our Galaxy and has super solar metal- 
licity (e.g., Peterson & Green 1998; Chaboyer et 
al. 1999; Stetson et al. 2003; Carney et al. 2005; 
King et al. 2005). This cluster, almost as old as 
the Galactic disk, is of paramount importance to 
study the time evolution of the disk properties. 
Apart from its age, NGC 6791 is interesting for 
his peculiar Horizontal Branch (HB), mostly com- 
posed by red stars, but with a (unusual) blue tail 
(Kaluzny & Udalski 1992; Liebert et al. 1994). Its 
study could be relevant for a number of issues, like 
e.g., the UV upturn in elliptical galaxies, the HB 
morphology and its connection with mass loss in 
globular clusters. 

Despite its peculiarities and its importance for 
Galactic formation studies, only one detailed work 
dealing with the chemical pattern in NGC 6791, 
based on modern fine abundance analysis and high 
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resolution spectroscopy, has been published so far: 
Peterson & Green (1998) analyzed the coolest and 
brightest (at V = 15.0, B — V = 0.48) blue HB 
star. Besides this case, really high S/N, high res- 
olution spectra can be obtained in acceptable ex- 
posure times only for the brightest, hence cooler 
giants in this cluster. At the very large metallicity 
of NGC 6791 these pose a great challenge to the 
observers (see next Section) and great care has to 
be taken to ensure the reliability of the analysis of 
their extremely crowded spectra. 

This paper is organized as follows: Sect. 2 
presents our data; atmospheric parameters and 
iron abundances from spectrum synthesis are de- 
scribed in Sect. 3; Sect. 4 provides the abundances 
of carbon and oxygen; a discussion and a summary 
are given in Sect. 5. 

2. Observational material 

The spectra of the giants of NGC 6791 are ex- 
tremely rich of lines, due to the rather cool tem- 
perature and the high metal content. To alleviate 
the analysis problems, we focused our attention on 
(fainter) stars on the red clump, which are warmer 
than the red giants. Even these spectra are actu- 
ally at least as rich in lines as that of the canonical 
very metal-rich giant /z Leo (see Figure 1). The 
long debate on the appropriate abundance to be 
attributed to this star (see e.g. Gratton & Sneden 
1990) emphasizes the difficulties of the derivation 
of correct abundances from the very line-rich spec- 
tra of cool giants. 

We chose our targets among the stars listed 
as cluster members by Friel et al. (1989) on the 
basis of their radial velocities (RVs). They are 
all confirmed members by our own measurements. 
All spectra were obtained with the high resolu- 
tion spectrograph SARG mounted at the Italian 
National Telescope Galileo (TNG) on Canary Is- 
lands. The resolution is R = 29, 000 and the wave- 
length coverage is 4620 - 7920 A. Individual spec- 
tra have been obtained (mostly in service mode) 
from November 2001 to September 2005, with ex- 
posures ranging from 1 to 1.5 hours each. For 
each star, the total exposure time ranges from 4.5 
to 8 hours. All spectra have been reduced using a 
standard IRAF 1 procedure for bias and flat field 

1 IRAF is distributed by the National Optical Astronomical 



correction, spectra extraction, and wavelength cal- 
ibration; the individual RVs were derived and the 
spectra were shifted to zero RV and averaged. Ta- 
ble 1 lists the identifications, coordinates and mag- 
nitudes for the stars, together with the heliocen- 
tric RVs obtained averaging all individual ones for 
each stars. Also shown is the S/N of the summed 
spectra measured near 6000 A by comparing the 
spectrum of the stars with those of \x Leo (see be- 
low) , that can be considered virtually noiseless for 
the purposes of this comparison. 

To cross-check our method of abundance deriva- 
tion, we also analyzed a spectrum of fi Leo with 
the same procedure. This star was selected be- 
cause it has stellar parameters and iron abun- 
dances similar to those expected in the NGC 6791 
sample. The spectrum of \x Leo was acquired us- 
ing the FEROS spectrograph at the ESO 1.5 m 
telescope at La Silla. The original spectrum has 
a resolution of R ~ 48, 000; however, in order to 
make the analysis as similar as possible to that 
of the stars of NGC 6791, this spectrum was de- 
graded at the same resolution of our SARG spec- 
tra. As evident from Figure 1, at this resolution 
the spectrum of fx Leo is indeed very similar to 
those of the stars of NGC 6791, except for the 
higher S/N, and very subtle differences in the line 
strengths. The fact that we are dealing with a 
different, lower resolution spectrum than the one 
analyzed by Gratton & Sneden (1990) and adopt a 
different solar model explain the slight differences 
in the results with that paper (see Sect. 3.4). 

3. Abundances in NGC6791 from synthe- 
sis of Fe lines 

In NGC 6791 even the spectra of red clump 
stars are so crowded that we deemed the tradi- 
tional analysis based on equivalent widths measure 
not entirely reliable, due to the large inherent un- 
certainties in both the location of the continuum 
level and the presence of blends (see Gratton & 
Sneden 1990 and Smith & Ruck 2000 for a simi- 
lar approach in the case of spectra of much higher 
resolution and S/N of \i Leo). 

We then derived iron abundances for our pro- 
gram stars by comparing the observed profiles for 



Observatory, which are operated by the Association of Uni- 
versities for Research in Astronomy, under contract with 
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a number of Fe lines to syntheses of small spectral 
regions (typical width ~ 2.5 A) around the cho- 
sen lines. This procedure allows to take fully into 
account the presence of blends. Furthermore, the 
correct positioning of the continuum level is much 
less a problem, since we may compare directly the 
highest points of the spectrum in the observed re- 
gions with those present in the synthetic spectra, 
insofar we trust the line lists used in our analysis. 
These lists were built taking appropriately into ac- 
count both the inclusion of all relevant lines and 
the quality of the gf values, carefully discussed in 
Gratton et al. (2003). 

These (extensive) line lists were obtained after 
careful comparisons with both the spectrum of the 
Sun and of HR 3627. This is a cool (~ 4200 K), 
very metal-rich ([Fe/H]~ +0.3) star. With this 
combination of parameters the strength of lines in 
HR 3627 is typically similar to (or even stronger 
than) that of the red clump stars of NGC 6791. 
We are fully confident that a line list well matching 
the spectrum of this star also provides sensible re- 
sults for our program stars. This lengthy prepara- 
tory work on HR 3627 was based on excellent ob- 
servational material: its spectrum, acquired with 
SARG, has both very high S/N (> 400) and reso- 
lution (R ~ 150, 000), so that even extremely faint 
possible contaminants could be detected and in- 
cluded in the line lists. We selected a number of Fe 
I and Fe II lines that were free from nearby strong 
features. We restricted to the spectral range from 
about 5500 to 7000 A, to avoid the low response of 
the spectrograph in the blue region and fringing or 
severe telluric contamination redward of the upper 
wavelength limit. Around each lines (within ±2 A) 
we extracted lines of neutral and singly ionized 
atomic species from the Kurucz database (Kurucz 
1995b). We also included molecular lines, in par- 
ticular of contaminant CN and hydrides. Lines 
unaccounted for in the Kurucz database and/or in 
the solar tables (Moore et al. 1966) were assumed 
to be Fe I lines with an excitation potential of 
3.5 eV. The transition probabilities of the Fe lines 
we wanted to synthetize were left untouched with 
respect to the line list we used in the EW analysis; 
the gf values of the nearby lines were adjusted one 
by one by matching the high resolution spectrum 
of HR 3627. For more details about HR 3627 and 
these line lists, see Carretta et al. (2004). 

The synthetic spectra were obtained using 



model atmospheres extracted from the grid of Ku- 
rucz (1995a). For consistency with other papers 
analyzing stars in open clusters (see e.g. Carretta 
et al. 2004), the models considered in this paper 
are those with the overshooting option included. 

The fitting of the synthetic spectra to the ob- 
served ones was done by eye. There is of course 
some arbitrariness in the eye fitting, since differ- 
ent weights can be attributed to the line cores or 
wings. However, we found that our eye fitting 
gives the same average abundances, but much less 
line-to-line scatter, than a fitting based on more 
objective criterions, such as a least square fitting 
to the data. The reason for the smaller scatter 
likely reflects a better estimate of the appropriate 
level of the local continuum, which is a free param- 
eter in the fitting, and it is a critical issue when 
determining abundances. 

The line parameters (oscillator strengths, damp- 
ing broadening) were obtained using the same pre- 
cepts adopted in Gratton et al. (2003). The same 
line parameters and microturbulent velocity were 
adopted for all stars. 

3.1. Atmospheric Parameters 

Effective temperatures (T e g) and surface grav- 
ities (log g) were obtained from the photometry, 
using B, V values from Montgomery et al. (1994) 2 
and K magnitudes from 2MASS (Cutri et al. 
2003). We derived the T eff 's from the V — K colors 
assuming the calibration by Alonso et al. (1999) 
and a reddening of E(B — V) = 0.15, in the mid- 
dle of literature determinations. We prefer not to 
use temperatures from the B — V colors because 
of their strong dependence on metal abundance. 
We have found a posteriori that the values of the 
effective temperatures corresponding to the B — V 
colors agree well with those from the V — K colors 
for the metal abundance of NGC 6791 determined 
in the present analysis. Also, we prefer not to 
use temperatures derived from line excitation be- 
cause for the lines considered here there is a quite 
strong degeneracy between effective temperatures 
and microturbulent velocities (see below). 

The calibration of the colour-temperature rela- 
tions of Alonso et al. (1999) was used for consis- 
tency with the analyses we are performing of sev- 

2 but other photometries are almost identical, see e.g. Stet- 
son et al. (2003) 
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eral other open clusters. It should be noticed that 
their calibration only extends up to [Fe/H]=+0.2, 
and application to the stars of NGC 6791 requires 
an extrapolation. However, the V — K color index 
used in this paper is only very marginally sensitive 
to metal abundance, so that errors in our tempera- 
tures cannot be large. Errors in our effective tem- 
peratures are mainly due to uncertainties in the 
assumed reddening: an error of AE(B — V) = 0.04 
magnitudes, a reasonable upper limit in the case 
of NGC 6791, implies an error of ~ 100 K in the 
assumed temperatures. For consistency, a similar 
approach was adopted for /z Leo, with the V — K 
color from Johnson (1966) photometry. Notice 
that this T e g for fi Leo is slightly cooler (by 50 K) 
than that derived by Gratton & Sneden (1990) us- 
ing the Infrared Flux Method. 

Surface gravities were obtained from the loca- 
tion of the stars in the color-magnitude diagram, 
assuming a distance modulus of (to— M)y = 13.45 
- in the middle of literature determinations - bolo- 
metric corrections from Alonso et al. (1999), and a 
mass of 0.9 M . Most of the errors in the surface 
gravities stem from uncertainties in the distance 
modulus: an error of (m — M)y = 0.5 mag, a 
reasonable upper limit, implies an error of 0.2 dex 
in the surface gravities. For fi Leo, we adopted 
the value of the gravity given by Gratton & Sne- 
den (1990) using a variety of methods (equilibrium 
of ionization for Fe, dissociation equilibrium for 
MgH, and pressure broadened lines): this value is 
in fact almost identical to that of the NGC 6791 
stars. 

Microturbulent velocities (w t ) were obtained by 
eliminating trends of abundances with respect to 
the expected line strength (Magain 1984) X = 
log gf — EP x <d exC7 where log gf is the oscilla- 
tor strength of the lines, EP the excitation poten- 
tial (in eV), and O exc = 5040/(0.86 x T off ) rep- 
resents the approximate temperature of the lay- 
ers where most of the lines form. This implies 
that abundances are roughly independent of line 
strength. We considered here the average values 
of the abundances derived for the individual lines 
of all the stars of NGC 6791, in order to reduce 
the scatter and better evidentiate possible trends. 
This means that the same value of the microturbu- 
lent velocity was adopted for all stars. The same 
value was also adopted in the analysis of /x Leo, 
and it is slightly smaller (by 0.15 km s _1 ) than 



what Gratton & Sneden (1990) and Smith & Ruck 
(2000) adopted. Uncertainties in these microtur- 
bulent velocities are approximately of 0.08 km s _1 
for a given effective temperature; this is obtained 
by modifying vt from its best value until the slope 
becomes equal to its statistical error. Note how- 
ever that there is a strong correlation between mi- 
croturbulent velocities and effective temperatures: 
adopting temperatures 100 K higher, we would 
have derived microturbulent velocities ~ 0.13 km 
s _1 larger. 

Final metallicities are obtained by interpolating 
in the Kurucz (1995a) grid of model atmospheres 
(with the overshooting option on) the model with 
the proper atmospheric parameters whose abun- 
dance matches that derived from Fe I lines. 

3.2. Abundances from individual lines 

Fig. 2 shows an example of the quality of the fit 
of three of the observed iron lines with synthetic 
spectra. The synthetic spectra have been com- 
puted with atmospheric parameters appropriate 
for the star, and Fe abundances of logn(Fe)=7.6, 
7.8, 8.0, 8.2, 8.4, 8.6, and 8.8. The lines consid- 
ered in the analysis are marked with a dash in 
Fig. 2. From these three comparisons we con- 
cluded for a best value of the Fe abundance of 
logn(Fe)=8.00. Note that other Fe I lines falling 
in the same spectral ranges of the program lines 
generally give Fe abundances in good agreement 
with those indicated by the lines selected in our 
analysis, although they were not actually used in 
the estimate of the best value for each star. 

Table 2 gives abundances for Fe lines 3 as ob- 
tained from the comparison of observations with 
spectral synthesis. Lines with expected line 
strength X > —5.6 were not used in the anal- 
ysis because there is a trend for these lines to give 
too low abundances. This is likely a reflection of 
the weight we gave in our abundance estimates to 
the cores of these lines. These form at very tiny 
optical depths, where the adopted model atmo- 
spheres are probably not adequate and deviations 
from LTE become important. Anyhow, such lines 
are strongly saturated, so that they would not be 

These abundances are by number; we use the usual spectro- 
scopic notations: logn(A) is the abundance (by number) 
of the element A in the scale where logn(H) = 12; [A/H] is 
the logarithmic ratio of the abundances of elements A and 
H in the star, minus the same quantity in the Sun. 
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good abundance indicators. Since the four target 
stars are very similar to each other, it is meaning- 
ful to average the results for the individual stars, 
and to derive the average abundance provided by 
each line. This is given in the last column of 
the table; the associated error bar is simply the 
dispersion of the mean. 

Fig. 3 shows graphically the lack of trends in 
this average Fe abundances with wavelength, EP, 
and line intensity parameter X, thus reinforcing 
the reliability of the derived abundances. 

Table 3 summarizes the abundances obtained 
for each star. The second, third and fourth 
columns give the values adopted for the atmo- 
spheric parameters (T e g, logg, [A/H]), while the 
microturbulent velocity is the same for all stars 
(v t =1.05 km s _1 ). In Columns 5 to 10 we 
give the average values of the abundances of Fe 
from neutral and singly ionized lines, along with 
the number of lines used in the analysis and the 
r.m.s. scatter of individual abundances. Finally, 
in the last two columns we give the [Fe/H] values. 
The reference solar abundances adopted here are 
logn(Fe) = 7.54 and 7.49 from neutral and singly 
ionized lines, respectively; those are the values we 
obtained from a solar analysis compatible with the 
present one (see Gratton et al. 2003). 

The last line of Table 3 gives the average abun- 
dance for the cluster. The average Fe abundance 
from neutral lines is [Fc/H] = +0.47 ±0.04 (r.m.s. 
of individual stars equal to 0.07 dex). The line-to- 
line scatter for each individual star is in the range 
0.11-0.18 dex. This is dominated by uncertainties 
still present in the location of the continuum level 
within the small spectral windows considered in 
this analysis. 

If we consider the average abundances for each 
line derived from the four spectra, the line-to-line 
scatter is 0.08 dex, as expected by assuming that 
the four results are independent of each other. For 
these average abundances, there is a small, not sig- 
nificant, trend with line excitation: A[Fe/H]/AEP 
= 0.009 ±0.012 dex/eV. This implies that temper- 
atures from line excitation are 38 ± 50 K higher 
than those adopted here. We deem this agreement 
as fully satisfying. 

There is also a small - but again not significant 
- difference between the abundances of Fe I and 
Fe II: [Fe/H]i - [Fe/H]n = 0.04 ± 0.07 dex. This 



small difference would have canceled, had we cho- 
sen temperatures about 50 K lower than adopted, 
or gravities about 0.15 dex larger. 

3.3. Errors on the derived abundances 

Table 4 gives the sensitivity of the derived 
abundances on the assumptions on the atmo- 
spheric parameters of the program stars (listed 
in Column 1). Column 2 gives the considered 
parameter variation, and Columns 3 and 4 the re- 
sulting variation in the abundances derived from 
Fe I and Fe II lines respectively. Column 5 lists 
the estimated value for the systematic error in our 
analysis for each of the parameters, and Columns 
6 and 7 the corresponding uncertainties in the 
abundances from Fe I and II lines, respectively. 

The total errors listed on the bottom line of 
Table 4 are obtained summing quadratically the 
contribution of each source of error, including the 
fitting error. This estimate of the total error is 
computed assuming a zero covariance between the 
effects of errors in the atmospheric parameters. 
In principle, this assumption is not strictly valid, 
since there are correlations between different pa- 
rameters. In practice, the Fe abundance is mainly 
a function of the adopted effective temperatures. 
In fact, adopting T c g values e.g. 200 K higher, we 
would have obtained higher microturbulent veloc- 
ities (by 0.26 km s _1 ), and Fe abundances smaller 
by about 0.15 dex. Notice that since the effect on 
microturbulent velocity is larger than the direct 
effect of temperature variations, the abundances 
get smaller with increasing model temperatures, at 
variance with the usual dependence for late type 
stars. Unluckily, given the correlation between ef- 
fective temperature and microturbulent velocity, 
we could not derive reliable effective temperatures 
from the excitation equilibrium. 

There are quite strong arguments favoring the 
solution adopted throughout this paper. In fact, 
had we adopted the hypothetical temperatures 
warmer by 200 K mentioned above (that would 
be obtained assuming a reddening of E(B — V) = 
0.23, a value near the upper limit in the analy- 
ses of NGC 6791), abundances from Fe n lines 
would have resulted much lower than those derived 
from Fe I lines, with an offset of 0.24 dex, rather 
than 0.04 dex. In principle, it should be possi- 
ble to compensate for such an offset by increasing 
the surface gravities by 0.5 dex (up to values of 
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logg ~ 2.8). However, this would require a dis- 
tance modulus of (m — M)v ~12.15, incompati- 
ble with the distances derived from the literature 
color-magnitude diagrams. The absolute magni- 
tude of core He-burning stars (i.e., the clump stars 
we are observing) is expected to be My ~ 1.2 
according to the models by Girardi & Salaris 
(2001). This corresponds to a distance modulus 
of (m — M)v ~ 13.4 — 13.5, in agreement with 
the value used in our analysis. Hence, an upper 
limit for the effective temperatures of the program 
stars in NGC 6791 is about 100 K higher than the 
adopted values (corresponding to a reddening of 
E(B — V) < 0.18, and a microturbulent velocity 
of Vt < 1.18 km s _1 ). This implies that the lower 
limit of the Fe abundance is [Fe/H]> +0.39 dex. 
Similar arguments can be used to define a robust 
upper limit of [Fe/H]< +0.55. 

Summarizing, the Fe abundance we derive for 
NGC 6791 is [Fc/H] = +0.47 + 0.04 + 0.08, where 
the first error value represents the random term, 
as derived from the star-to-star scatter, and the 
second error value the systematic errors, due to 
the assumptions on reddening, distance modulus, 
mass, temperature scales, etc. 

3.4. A comparison with /j, Leo 

To further assess the soundness of the results 
of our spectrum synthesis analysis, we also ana- 
lyzed the well known, metal-rich field star \x Leo, 
strictly using the same criteria. The average Fe 
abundance we derived is logn(Fe) = 7.92 ± 0.03 
(29 lines, r.m.s.=0.15 dex) from neutral lines, and 
logn(Fe) = 7.79 ± 0.04 (6 lines, r.m.s.=0.10 dex) 
from singly ionized lines. These values correspond 
to [Fe/H]=+0.38±0.03 and [Fe/H]=+0.30 ± 0.04 
respectively. The difference in the abundances 
provided by neutral and singly ionized lines is only 
marginally larger than the internal error, and is 
similar to the corresponding difference found for 
the stars in NGC 6791. 

This value of the Fe abundance of fi Leo can be 
compared with literature estimates. We will con- 
sider only abundances based on very high quality 
spectra (R ~ 100, 000, S/N > 200) and compar- 
isons with synthetic spectra, since use of equiva- 
lent widths is not reliable for such line-rich spec- 
tra. Gratton & Sneden (1990) derived logn(Fe) = 
7.97 ± 0.02 ± 0.15 (36 lines, r.m.s.=0.12 dex). 
Within the internal errors, the present estimate 



of the abundance agrees with that of Gratton & 
Sneden (1990). The [Fe/H] value given by Gratton 
& Sneden (1990) is slightly lower than that esti- 
mated in this paper ([Fe/H]=+0.34 ± 0.03), due 
to the use of a different solar abundance, derived 
using the Holweger & Muller (1974) model atmo- 
sphere, rather than that obtained from the same 
grid considered for the red giants. More recently, 
Smith & Ruck (2000) used a very similar technique 
and set of atmospheric parameters, and obtained 
a lower value of logn(Fe) = 7.79 ± 0.03 (internal 
error). This value is about 0.13 dex lower than 
the present one and the difference can be ascribed 
to the larger value of the microturbulent velocity 
adopted by Smith & Ruck (2000) (1.22 km s" 1 
rather than 1.05 km s _1 : see Table 4). These 
comparisons show that systematic errors in our 
analysis are within the adopted error bar. 

By comparing the abundances in NGC 6791 
with those of \i Leo, we conclude that the for- 
mer has an Fe abundance larger than the latter 
by 0.09 ± 0.05, where the error here is simply 
the quadratic sum of the internal errors, since the 
other uncertainties should cancel out in the dif- 
ferential abundance analysis. This small positive 
difference well agrees with the eye impression from 
Figure 1. 

4. Abundances of oxygen and carbon 

We estimated the abundance of oxygen from 
the [O i] line at 6300.31 A. The second, weaker, 
[O i] line at 6363.79 A is barely measurable in our 
spectra, and the three O lines near 7775 A (that 
are not the best suited for this kind of stars) fall 
in a region where strong interference fringes make 
it difficult to extract accurate information from 
spectra. 

The 6300.31 A line often happens to be strongly 
affected by telluric contamination. However, in 
the case of NGC 6791, the cluster RV relative to 
Earth motion makes the [O i] line non contami- 
nated by telluric absorptions in the majority of our 
spectra. Compromising between a slightly lower 
S/N and a lower data manipulation, we decided 
to average only the undisturbed spectra for each 
star, and measure the O abundance from them. 

In order to derive reliable O abundances, a care- 
ful synthesis of the [O i] lines is necessary, includ- 
ing not only the contribution of the nearby Ni line 
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but also the relevant coupling with C abundances, 
and the contamination by CN lines. While a full 
discussion of CNO abundances is deferred to a 
forthcoming paper, here we preliminarily estimate 
the C abundances from the spectral synthesis of 
the C2 molecular features at 5086 A. Fig. 4 shows 
an example of matching of synthetic spectra for 
star 3019 in the forbidden [O 1] line and in the 
C2 regions. The synthetic spectra were computed 
assuming [N/Fe]=-0.2, and 12 C/ 13 C=8. The lat- 
ter value is adequate for low mass evolved giants, 
while the former one provides CN features in rea- 
sonable agreement with observations for the best 
C abundance. 

Combining the C and O abundances provided 
by these two abundance indicators, we find aver- 
age abundance ratios of [C/Fe] = —0.23 dex (a 
rather normal value for clump stars), and [O/Fe] 
= —0.32 dex. The comparisons with synthetic 
spectra clearly show that the O abundance can- 
not be much larger than this value. It cannot be 
much lower either, otherwise the features due to 
C-bearing molecules would be much stronger (the 
spectra of the stars clearly shows that 0^>C). 

As a comparison, Gratton & Snedcn (1990) de- 
rived for fi Leo [C/Fc] = [O/Fe] = -0.15 dex. 
NGC 6791 is then slightly more deficient in both 
O and C with respect to /x Leo, along the stan- 
dard trend of more pronounced deficiency of these 
elements with increasing metal abundance (see 
Bensby et al. 2005; Andersson & Edvardsson 1994) 

5. Discussion and summary 

The first detailed study of NGC 6791 of which 
we are aware of is by Kinman (1965), who pub- 
lished a photographic color-magnitude diagram 
and derived information on membership and in- 
trinsic colors (i.e., reddening) from low resolution 
spectra. 

There is a general agreement that this cluster 
is at about the same Galactocentric distance as 
the Sun, is about twice as old, and about twice 
as metal-rich. Nevertheless, there are significant 
differences between properties derived by different 
authors: (m - M) ~ 12.6 - 13.6, E(B - V) ~ 
0.09-0.23, age - 7-12 Gyr, [Fc/H] - 0.1-0.4 (for 
a recent review see Stetson et al. 2003). NGC 6791 
has also been observed with the ACS on board the 
Hubble Space Telescope (King et al. 2005; Bcdin 



et al. 2005), reaching almost the hydrogen-burning 
limit on the main sequence and defining the White 
Dwarfs cooling sequence. 

Spectroscopic studies of NGC 6791 are fewer, 
but present interesting results. Since the pioneer- 
ing work by Spinrad & Taylor (1971) the cluster 
has been recognized to be more metal-rich than 
the Sun; they found [M/H] = +0.75, but with a 
method that attributes [M/H] ~ +0.6 to NGC 188 
and M 67, both presently known to have nearly so- 
lar metallicity. 

More recently, the red HB (that we call red 
clump) has been studied by low resolution spec- 
troscopy by Hufnagel et al. (1995); they wanted to 
investigate possible correlations between CH and 
CN band strengths and found none, at variance 
with similar studies on globular clusters. 

Worthey & Jowett (2003) took low resolution 
spectra of 23 K giants in NGC 6791 and derived, 
using Lick/IDS indices, [Fe/H] = +0.320 + 0.023 + 
<r S ys, where a sys , the systematic error attached to 
their metallicity scale, is unknown. We have only 
one star in common (3019, their R25) for which 
the quoted abundance is [Fe/H] = +0.341. 

All the four stars analyzed here were observed 
at low resolution by Friel et al. (1989), where 
we adopted identification and membership sta- 
tus from. The same group (Friel & Janes 1993) 
later published their metallicities: [Fe/H] values 
are +0.32 ± 0.28, +0.23 ± 0.33, +0.39 ± 0.15, 
+0.41 + 0.21 for stars 2014, 3009, 3019, and SE49, 
respectively. These values (that would give an 
average metallicity of about +0.34) do not dif- 
fer much from our new ones, considering that the 
scale used by Friel & Janes (1993) gives on average 
metal-poorer values than others. The final result 
of their analysis depends on the calibration of in- 
dices using standard stars with metallicity derived 
by several different sources. Due to a recent revi- 
sion of their abundance scale (Friel et al. 2002), 
these abundances were lowered to +0.12 ± 0.13, 
+0.11 + 0.13, +0.15 + 0.11 and +0.16 + 0.11 for the 
same stars; for their total sample of 39 stars they 
derived an average metallicity of +0.11, <r = 0.10 
dex. 

To date, the only published analysis based on 
high resolution spectroscopy remains that by Pe- 
terson & Green (1998). They obtained 4 hours 
of integration at the 4 m Mayall telescope at Kitt 
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Peak on a blue HB star that is a cluster member 
both by proper motion and RV. The summed spec- 
trum has S/N and resolution lower than ours (S/N 
~ 30 per pixel, R = 20000), but is less affected by 
blends because of the much higher temperature 
(about 7300 K). They employed both EWs (for Fe 
and atmospheric parameters determination) and 
spectrum synthesis (for all the other measured ele- 
ments). They discussed the use of different scales, 
e.g. for logg/'s and solar reference abundances. 
On the scale they chose to adopt, their star has 
[Fe/H] = +0.37 (±0.10) dcx; since their solar Fe 
is logn = 7.67 and ours is 7.54, their value corre- 
sponds to [Fe/H] = +0.50 on our scale. 

Taken at face value, the agreement with our re- 
sults appears excellent, but it doesn't take into ac- 
count many differences like e.g., logg/'s or model 
atmospheres and the like. For instance, we have 4 
Fe lines in common with Peterson & Green (1998), 
1 of Fe I and 3 of Fe n (see Table 2 and their 
table 1), and the average difference in log<?/'s 
is about 0.29 (our values minus theirs). They 
also measured the abundances of several elements 
by means of extensive spectrum synthesis, finding 
strong overabundances with respect to the solar 
ratio for N and Na (+0.5 or +0.4) and slightly 
lower ones for Mg and Si (+0.2) or Eu (+0.1). All 
the other elements, including C and O, display 
solar-scaled abundances. In our analysis we have 
measured so far only C and O abundances, finding 
for them sub-solar ratios. A complete comparison 
between the two analyses is beyond the purpose 
of the present paper, also because the chosen tar- 
gets are quite different in properties, but it will be 
done once we have measured all other elements. 

Summarizing, we have derived the metallicity 
of NGC 6791 analyzing the high-resolution, high 
S/N spectra of four red clump stars, finding an 
average [Fe/H] = +0.47 ± 0.04 ± 0.08, where the 
two error values represent the random term and 
the systematic one, respectively. 

This result has been obtained using spectrum 
synthesis of Fe I and Fe n lines, and adopting at- 
mospheric parameters from the photometry. The 
soundness of our analysis has also been checked 
relatively to a well studied, metal-rich field giant 
star, /i Leo. 

We have measured O and C abundances using 
synthesis of the [O i] line at 6300.31 A and of the 
C2 molecular feature at 5086 A, finding [O/Fe] = 



—0.32 and [C/Fe] = —0.23 dex, on average, for 
NGC 6791. 

NGC 6791 represents a challenge to any model 
for the formation of the Galactic disk, a possible 
missing link between globular and open clusters, 
a genuine puzzle that modern abundance analysis 
based on high resolution spectroscopy will help to 
unveil. 

We thank the TNG staff who helped collecting 
the spectra. We gratefully acknowledge the use 
of the BDA, created and updated by J.-C. Mer- 
milliod for many years. This project has been 
partially supported by the Italian MIUR, under 
PRIN 2003029437. 
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Table 1 

Information on the four observed stars. 



ID 


ID 


RA 


Dec 


V 


B-V 


K 


RV 


S/N 


(FLJ) 


(BDA) 


(2000) 


(2000) 






2MASS 


km s -1 




2014 


2562 


19 21 01.1 


37 46 39.6 


14.563 


1.403 


11.468 


-48.60 


60 


3009 


3363 


19 20 56.5 


37 44 33.7 


14.638 


1.387 


11.556 


-48.26 


50 


3019 


3328 


19 20 56.2 


37 43 07.7 


14.675 


1.334 


11.586 


-46.10 


85 


SE49 


3926 


19 20 53.1 


37 45 33.4 


14.540 


1.346 


11.513 


-45.63 


40 



References. — FLJ: Friel et al. (1989); BDA: Mermilliod (1995); B, B-V: Montgomery 
et al. (1994); K: Cutri ct al. (2003); RV: this paper 
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Table 2 

Abundances from individual Fe lines 



Wavel. 

(A) 


EP 

(eV) 


loe q f 


/i Leo 


2014 


3009 


3019 


SE49 


Mean 










Fe I 










5560.22 


4.43 


-1.10 


7.71 




8.00 


8.06 


7.91 


7.99 ±0.04 


5577.03 


5.03 


-1.49 


7.76 


7.91 


8.24 


7.93 


7.86 


7.98 ±0.08 


5618.64 


4.21 


-1.34 


7.67 


7.74 


8.01 


8.02 


8.07 


7.96 ±0.08 


5619.61 


4.39 


-1.49 


7.70 


8.09 


8.19 


7.87 


8.04 


8.05 ±0.06 


5650.00 


5.10 


-0.80 


8.04 


8.12 


8.17 


8.27 


7.81 


8.09 ±0.10 


5651.48 


4.47 


-1.79 


7.98 


8.03 


8.37 


7.88 


8.01 


8.07 ±0.10 


5661.35 


4.28 


-1.83 


8.04 


8.03 


7.83 


8.13 


8.03 


8.01 ±0.06 


5784.67 


3.40 


-2.53 


7.84 


7.95 


7.91 


7.88 


7.87 


7.90 ±0.02 


6078.50 


4.79 


-0.38 


8.02 


7.74 


8.21 


7.85 


8.17 


7.99 ±0.12 


6079.02 


4.65 


-0.97 


7.86 


7.99 


8.19 


8.07 


8.15 


8.10 ±0.04 


6082.72 


2.22 


-3.53 


7.89 


7.91 


7.84 


7.87 


7.97 


7.90 ± 0.03 


6089.57 


5.02 


-0.87 


7.64 


7.93 


7.97 


7.87 


7.62 


7.85 ±0.08 


6093.65 


4.61 


-1.32 


7.74 




8.08 


7.97 


7.83 


7.96 ± 0.06 


6096.67 


3.98 


-1.76 


7.86 


8.04 


8.00 


7.97 


7.81 


7.96 ±0.05 


6098.25 


4.56 


-1.81 


7.82 


8.12 


8.16 


8.08 


7.80 


8.04 ± 0.08 


6137.00 


2.20 


-2.91 


7.92 


7.96 


8.25 


8.04 


7.95 


8.05 ±0.07 


6270.23 


2.86 


-2.55 


8.09 


8.01 


8.04 


8.16 


8.07 


8.07 ±0.04 


6380.75 


4.19 


-1.34 


7.92 




8.16 


8.06 


7.87 


8.03 ± 0.08 


6481.88 


2.28 


—2.94 


8.11 


7.86 


8.25 


8.00 


7.94 


8.01 ± 0.08 


6498.94 


0.96 


-4.66 


7.91 


7.87 


8.07 


7.86 


7.91 


7.93 ±0.05 


6518.37 


2.83 


-2.56 


8.24 


8.16 


8.14 


8.06 


8.47 


8.21 ±0.09 


6533.94 


4.56 


-1.28 


8.10 


7.83 


8.39 




8.24 


8.15 ±0.14 


6574.25 


0.99 


-4.96 


8.10 


7.98 


8.11 


7.98 


7.94 


8.00 ± 0.04 


6581.22 


1.48 


-4.68 


7.98 


7.77 


8.05 


7.89 


7.97 


7.92 ± 0.06 


6609.12 


2.56 


-2.65 


8.11 


7.96 


7.95 


8.05 


8.24 


8.05 ±0.06 


6633.76 


4.56 


-0.81 


8.07 


7.94 


8.00 


7.86 


8.27 


8.02 ±0.09 


6713.75 


4.79 


-1.41 


7.93 


7.77 


8.01 


7.98 


8.20 


7.99 ± 0.08 


6725.36 


4.10 


-2.21 


7.83 


7.73 


8.18 


7.98 


8.22 


8.03 ±0.11 


6726.67 


4.61 


-1.05 


7.90 


7.94 


8.19 


8.02 


8.05 


8.05 ±0.06 










Fe II 










6247.56 


3.87 


-2.32 


7.74 


8.01 


7.55 




8.50 


8.02 ±0.24 


6369.46 


3.89 


-4.21 


7.80 


8.02 


7.75 


7.93 


7.69 


7.85 ±0.12 


6383.72 


5.55 


-2.09 


7.94 


7.88 


8.16 


7.85 


8.24 


8.03 ±0.10 


6416.93 


3.89 


-2.70 


7.67 


7.94 


8.02 


8.05 


7.70 


7.93 ± 0.08 


6432.68 


2.89 


-3.58 


7.84 


7.73 


7.67 


7.62 


8.32 


7.84 ±0.16 


6456.39 


3.90 


-2.10 


7.75 


7.81 


8.07 


7.55 


8.12 


7.89 ±0.13 
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Table 3 

Summary of atmospheric parameters and abundance derivations; v t is 1.05 km s" 1 for all 
stars. Last line gives the cluster mean iron abundance. 



Star 


Teff 




[A/H] 




log n(Fe 


i) 




log n(Fe 


n) 


[Fe/H]l 


[Fe/H]n 




(K) 






n 


mean 


r.m.s. 


n 


mean 


r.m.s. 






fi Leo 


4490 


2.30 


0.47 


29 


7.92 


0.15 


6 


7.79 


0.10 


+0.38 


+0.30 


2014 


4463 


2.30 


0.50 


26 


7.94 


0.12 


6 


7.90 


0.11 


+0.40 


+0.41 


3009 


4473 


2.33 


0.50 


29 


8.10 


0.14 


6 


7.87 


0.24 


+0.56 


+0.38 


3019 


4468 


2.35 


0.50 


28 


7.99 


0.11 


5 


7.80 


0.21 


+0.45 


+0.36 


SE49 


4512 


2.32 


0.50 


29 


8.01 


0.18 


6 


8.09 


0.33 


+0.47 


+0.60 


(cluster) 










8.01 


0.07 




7.92 


0.08 


+0.47 


+0.43 
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Table 4 

Sensitivity and errors due to uncertainties in the atmospheric parameters 



Parameter 


Change 


A[Fe/Hi] 


A[Fe/H]n 


Error 


[Fe/H]i 


[Fe/H]n 


T c ff 


+100 K 


+0.016 


-0.158 


100 K 


+0.016 


-0.158 


logg 


+0.3 dex 


+0.036 


+0.184 


0.2 dex 


+0.024 


+0.123 


[A/H] 


+0.1 dex 


+0.024 


+0.044 


0.06 dex 


+0.014 


+0.026 


Vt 


+0.2 km s" 1 


-0.124 


-0.047 


0.13 km s- 1 


-0.075 


-0.028 


Fitting error 








0.16 dex 


+0.026 


+0.099 


total 










0.086 


0.226 



Note. — Fitting error is the line-to-line uncertainty related to the procedure of attributing 
an abundance to each individual line. It includes random errors due to the photometric 
uncertainties on the line profile and to the estimate of the local continuum level. 
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NGC6791 ■ SE49 




Q I I I I I I I I I l I I I I I I I 

6080 6082 6084 6086 6088 6090 

Wavelength (A) 

Fig. 1. — Small portion of the spectra of our target stars compared with the one of \i Leo, degraded to the 
same resolution. 
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Fig. 2. — Comparisons between the observed spectrum for a star in NGC 6791 (star 3009: thick line) and 
synthetic spectra (thin lines), for a few of the spectral regions containing Fe lines considered in the present 
analysis. 
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Fig. 3. — Runs of the abundances obtained by averaging the results for each individual Fe I line over the 
four program stars, with wavelength, excitation potential EP 1 and the expected line intensity parameter 
X = log gf — EP * 5040/(0.86 * T c ff). The error bars for the individual lines are the maxima between the 
spread of the mean over the four stars, and the r.m.s. of this quantity (0.08 dex). The average value of the 
abundance for NGC 6791 is indicated by the dashed horizontal line, while the formal best fit line based on 
Fe I is indicated by the solid line. There is no significant trend in any of the panels of this figure. 
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Fig. 4. — Comparison between the observed spectrum for star 3019 in NGC 6791 (thick line) and synthetic 
spectra (thin lines) for the [O i] line at 6300 A (upper panel) and the C2 molecular feature at 5086 A(lower 
panel). The synthetic spectra were computed assuming [N/Fe]=— 0.2, and 12 C/ 13 C=8. In the upper panel 
we adopted [C/Fe]=-0.25, and [O/Fe]=-0.65, -0.55, -0.45, -0.35, -0.25, -0.15, -0.05. In the lower panel 
we adopted [O/Fe]=-0.35, and [C/Fe]=-0.40, -0.35, -0.30, -0.25, -0.20, -0.15, -0.10. 
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